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Express MaU No. EV3366i6999US 

METHOD AND SYSTEM FOR DETERMINATION OF MOLECULAR 
INTERACTION PARAMETERS 

5 BACKGROUND OF THE INVENTION 

Field of the invention 

The present invention relates to the determination of kinetic parameters 
for molecular interactions, and more particularly to a method for determining 
kinetic parameters for the interaction between a molecule immobilised to a sdld 
10 support surface and a binding partner to the molecule in solution. The invention 
also relates to an anaiytical system, a computer program product and a 
computer .system for performing the method. 

Description of the related art 

15 Analytical sensor systems that can monitor Interactions between 

molecules, such as biomolecules, in real time are gaining increasing Interest. 
. These sehsor systems, usually referred to as interacjion analysis sensor 
systems or blospecific interaction analysis sensor syst^ras-. are often based on 
optical biosensors and affinity analysis and offer a rgpid way to detfermine in 

.20 real time inter alia equilibrium and rate constants without the need to label the 
interacting molecules. They have been used In the study of a variety of 
biomolecules, Including proteins, nucleic acids, lipids and carisohydrates. In 
these systems, a sensor surface having one of the molecular reactants 
immobilised thereto is contacted with a solution containing the other reactant. 

25 either by providing a flow of the solution past the sensor surface, or In a cuvette 
or the like, and binding interactions at the surface are detected. 

Conventionally, to determine, for example, association and dissociation 
rate constants (kg and k^, respectively) for the interaction between two 
Interacting molecules, one of the molecules, often refen-ed to as the ligand, is 

30 immobilised to a sensor surface and the other molecule, often referred to as the 
analyte. is provided in solution at several different known concentrations. Each 
concentration, or sample, of the analyte is then contacted with the sensor 
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surface, either in a laminar flow past the sensor surface or In a cuvette or the ' 
like to pemrjit association of the analyte to the sensor surface. After a' sample 
has been brought to contact the sensor surface, the surface is contacted with a 
solution free from analyte, usually buffer, to permit dissociation of the analyte 
from the immobilised iigand. During these association and dissociation phases, 
the amount of binding of analyte to the surface Is continuously detected and the 
binding data is collected. Before contacting the sensor surface with sample of a 
new analyte concentration, the Iigand surface Is restored or "regenerated" by 
treating the surface with a regeneration solution capable of removing any bound 
analyte while not destroying the Iigand. In that way, all the different samples will 
contact essentially one and the same iigand surface as far as Iigand density is 
concerned, the association and dissociation rate constants can then be • 
obtained from the collected binding data by fitting the data to mathematical 
descriptions of interaction models In the form of differential equations. Usually, 
the binding.data for all the samples are used in the same fit, a procedure ' 
refen-ed to as global fitting. From the determined association and dissociation 
rate constants kg and k^j. the equilibrium constant, Kq. and the affinity constant 
Ka (KA=1/KD)-of the interaction can in turn 66 calculated. Alternatively, 
provided that the interaction reaches steady state during the association phase, 
the equilibrium constant can be directly obtained from the binding data without 
fitting. 

Problems may arise, however, when the ligands are covalently 
immobilised to the sensor surface and suitable regeneration conditions are 
difficult to find. Renewed binding of the Iigand via an Immobiilsed capture agent 
before the contacting with each new analyte concentration could then be an 
alternative, but has the disadvantage of consuming large amounts of Iigand for 
the detennination. 

Detennination of equilibrium constants by a titration procedure without 
• the requirement for regeneration of the immobilised ligands is described by 
Schuck. P.. et al. (1998) Anal. Biochem. 265, 79-91. The sample is continuously 
circulated In a closed loop over two sensor spots of a commercial surface 
plasmon resonance biosensor. One of the sensor spots is functionalized wiUi an 



immobilised iigand for a soluble analyte in the sample, and the other sensor 
spot serves as a reference surface. A binding isotherm for the Interacting 
molecules is obtained by stepwise equilibrium titration of the analyte into the 
circulating loop, i.e. the sensor spots are sequentially contacted with stepwise 
5 Increased concentrations of analyte until equilibrium is attained for each 
concentration. This equilibrium titration is said to be especially useful for the. 
determination of binding constants in high-affinity systems since it eliminates 
the need for interpretation of binding idnetics and thereby problems that may 
arise from mass transport limitations. 

10 A similar stepwise equilibrium titration procedure is described for a 

cuvette-based biosensor design in Hall. D. R., and Winzor. D. J. .(1997) Anal. 
Biochem. 244. ^52-^e0. 

Also Myszka, D. G.". et al. (1998) Anal. Biochem. 265. 326-330 discloses 
equilibrium analysis of high affinity interactions using a surface plasmon 

15 resonance-based biosensor. In this approach, the time available to collect 
association phase data is increased by placing the analyte directly into the 
running buffer. Complete equilibrium binding profiles were generated without a 
regeneration step by changing the concentration of analyte and allowing the 
surface reactions to reequilibrate. Analyte concentrations were also decreased 

20 to demonstrate that the binding reactions were fully reversible. In this way. 
equilibrium dissociation constants for very high affinity interactions could.be 
detennlned. 

Shank-Retzlaff, M. L.. and Sligar. S. G. (2000) Anal. Chem. 72, 421 2- 
4220 describes a one-step method for detemiining kinetic rates and equilibrium 

25 binding affinities by a technique called analyte-gradient surface plasmon 

resonance SPR (AG-SPR) which eliminates the need for surface regeneration. 
A gradient of analyte is passed over the sensor surface under continuous-flow 
conditions so that the concentration of analyte increases llneariy with time. The 
rate at which analyte binds to Immobilised llgands on the sensor surface is 

30 monitored by monitoring the change in the surface plasmon resonance as the 
analyte concentration Increases. Kinetic rates are determined by fitting the data 
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to a two-compartment model for the molecular interaction which permits use 
also for systems under mass transport limitations. 

Titration procedures to detemiine binding capacity and regeneration 
conditions have been proposed by Karisson, R.. et al.. Biasymposium 1998, 
Edinburgh 2-4 September 1998, 5th Biasymposium in Japan, Tokyo, November 
5-6, 1998. 

US-A1-2003/014365 discloses a method for determining interaction 
parameters, including rate constants, between an analyte and a ligand 
immobilised to a sensor surface, where a measurement can be perfomied 
several times in succession, e.g. In a cuvette, with a stepwise modification of 
the analyte concentration each time. The measurements need not be completed 
to equilibrium t)ut can be ihtenupted earlier and the concentration of the analyte 
can be raised or lowered. Separate fitting of the binding curve part 
coreespondlng to each analyte concentration is used to detemiine- respective 
Initial binding rates, from which the association and dissociation fate constants 
for the interaction are then determined. 

Thus, while titration procedures in combination with biosensors and 
affinity analysis are Itnown perse in the art for determining equilibrium binding 
affinities, the use of multiple titrations to determine Idnetic rate constants seems 
fo be disclosed only in the above-mentioned publications Shank-Retzlaff, M. L, 
and Sllgar. S. G. (2000). and US-A1 -2003/014366. The method according to 
Shank-Retzlaff. M. L... and Sllgar. S. G. suffers from the limitation that it requites 
the use of a continuous gradient ofihe analyte. The method of US-A1- 
2003/014365. on the other hand, is disadvantageous jn that initial binding rates 
are frequently lower than the kinetic binding rate due to transport limitations, 
making initial binding rates unreliable for kinetic analysis. Further, this type of 
evaluation restricts the injection order of analyte concentrations. 

From the prior art It may therefore be concluded that for detennining 
kinetic rates for molecular interactions using systems based on biosensors and 
affinity analysis, it is necessary to regenerate the Immobilised ligand prior to 
contacting the sensor surface with a different concentration of analyte to 
thereby present essentially one and the same ligand surface to each analyte 



concentration, unless (i) a continuous gradient of tlie analyte is used, or (ii) 
initial binding rates are detennined In systems free firom mass transport 
limitations. 

It Is an object of the present invention to provide a sensor-based method 
5 for detemiining chemical interaction parameters, including l<inetic rate 
constants, by stepwise titration, which method obviates regeneration 
procedures while permitting measurements under mass transport limitation. 

BRIEF SUMMARY OF THE INVENTION 
10 According to the present invention, it has surprisingly been found that the 

above- and other objects and advantages can achieved by an affinity sensor- 
based method for detennining chemical interaction parameters, which method 
comprises perfomning, in one and the same experimental cycle, stepwise 
changes of the analyte concentration without intermediate regeneration of the 

•15 sensor surface, and determining the interaction parameters using global 

analysis of the whole set of detected binding data with regard, to a kinetic model 
for the molecular Interaction that allows for mass transport limitation, such as 
e.g. the afore-mentioned two-compartment model. In this method, Idnetlc data 
can thus be obtained in spite pf mass transport limitations ahd the different 

20 analyte concentrations may be used in any order. * * 

Therefore, In one aspect, the present invention provides a method of 
detennining Idnetic parameters for a- reversible molecular Interaction between a 
ligand Immobilised to a sensor surface and a bihding partner to the llgand In 
solution, comprising the steps ofr 

25 a) sequentially, without intermediate regeneration or renewal of the 

immobilised ligand, flowing a plurality of fluid volumes containing different 
l<nown concentrations of the binding partner over the solid support surface to 
pemiit association of binding partner to the Immobilised ligand, 

b) flowing over the solid support surface a fluid volume free from binding 

30 partner to pemiit dissociation of binding partner from the ligand, 
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c) monitoring during steps a) and b) the momentary amount of binding 
. partner bound to the solid support surface related to time and solution 

concentration of binding partner arid collecting the binding data, and 

d) determining the kinetic parameters by globally fitting a predetenmined 
5 kinetic model for the interaction between the binding partner and the 

immobilised ligand to the collected binding data, which model allows for mass 
transport limitation at the solid support surface. 

The solid support surface is preferably a sensor surface, i.e. one which 
produces a detectable signal In response to a binding interaction at the surface. 
10 The temi "monitoring" as used herein means that detection is perfonned 

at at least a plurality of times during steps a) and b), preferably at a large 
number of times. 

In another aspect, the present invention provides a method for 
detemriining the concentrations of binding partner in a plurality of binding 
15 partner-containing samples. 

In still another.aspect. the present invention provides a method for 
detemriining kinetic parameters for a plurality of different binding partners. 

In another aspect, the present Invention provides.an analytical System for 
• studying molecular interactions, which system comprises data processing 
20 means for performing at least one of the methods. 

In yet another aspect, the present invention provides a computer 
program comprising program code means for perfonning at least one of the 
methods. . 

In still another aspect, the present invention provides a computer 
25 program product comprising program code means stored on a computer 
readable medium or can-led on an electrical or optical signal for perfomiirig at 
least one of the methods. 

In another aspect, the present invention provides a computer system 
containing a computer program comprising program code means for perfonning 
30 at least one of tiie methods. 
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other advantages, novel features and objects of the invention will 
become apparent from the following detailed description of the invention when 
considered In conjunction with the accompanying dravirings. 



5 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 1s a schematic side view of a biosensor system based on 
surface plasmon resonance (SPR). 

Figure 2 is a representative sensorgram with a binding curve having 
association and dissociation phases. 
10 Figure 3 is a representative sensorgram obtained by sequential injections 

of a kinase inhibitor (staurosporin) over a surface with immobilised kinase 
.according to the method of the present invention. 

Figure 4 is the sensorgram of Figure 3 having overlaid thereon the • 
binding curve resulting from global fitting. 
15 Figure 5 is an overlay plot of obsen/ed and fitted binding data for 

sequential injections of a carbonic anhydrase Inhibitor (acetazolemide) over a 
surface with immobilised carbonic anhydrase. 

Figure 6 is a similar overlay plot to that of Figure 5 for sequential 
• • injections of another carbonic anhydrase Inhibitor (azosulfamide). ■ 
20 Figure 7 is a similar overlay plot to those of Figures 5 and 6 for 

sequential injections of still another carbonic anhydrase inhibitor 
(benzenesulfonamide). 

Figure 8 is an overlay plot of observed and fitted binding data for 
sequential injections of a thrombin inhibitor (melagatran) over a surface with 
25 immobilised thrombin. 

Figure 9 is an overlay plot of four binding curves and fitted data, three • 
curves (to the left) representing sensorgrams with single injections of analyte, 
and one curve (to the right) representing a sensorgram for sequential injections 
of analyte. 
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DETAILED DESCRIPTION OF THE INVENTION 

Unless defined otheronse, all technical and scientific temis used herein 
have the same meaning as commonly understood by a person skilled in the art 
related to this Invention. Also, the singular forms "a", "an", and "the" are meant 
5 to include plural reference unless it Is stated othennrise. 

As mentioned above, the present invention relates to the determination 
of molecular interaction parameters, including kinetic rate cohstants, for the 
interaction between a molecule immobilised to a solid support surface and a 
binding partner to the molecule In solution by a novel titration type method. 

10 preferably in combination with sensor based technology to study the molecular 
interactions and present the results in real time, as the interactions progress.- • 
Before describing the present Invention in more detail, however, the general 
context in which the invention is used will be described. ' • 

Chemical sensors or biosensors are typically based on label-free 

15 techniques, detecting a change in a property of a sensor surface, such as e.g. 
mass, refractive index, or thickness for the immobilised layer, but there are also 
sensors relying on some kind of labelling. Typical sensor detection techniques 
include, but are not limited to, mass detection methods, such as piezoelectric, 
optical, thermo-optical and surface acoustic wave (SAW) device methods, and 

20 electrochemical methods, such as potentiometric, conductometric, 

amperometric and capacitance/impedance methods. With regard to optical 
detection methods, representative methods include those that detect mass 
surface concentration, such as reflection-optical methods, including both 
external and internal reflection mefliods, angle, wavelength, polarization, or 

25 phase resolved, for example evanescent wave ellipsometry and evanescent 
wave spectroscopy (EWS, or Internal Reflection Spectroscopy), both may 
include evanescent field enhancement via surface plasmon resonance (SPR). 
Brewster angle refiractometry, critical angle refractometry, frustrated total 
reflection (FTR), scattered total internal reflection (STIR) - which may include 

30 scatter enhancing labels, optical wave guide sensors: external reflection 
imaging, evanescent wave-based imaging such as critical angle resolved 
imaging, Brewster angle resolved Imaging, SPR-angle resolved imaging, and 
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the like. Further, photometric and imaging/mlcrx)scopy methods, "per se" or 
combined with reflection methods, based on for example surface enhanced 
Raman spectroscopy (SERS), surface enhanced resonance Raman 
spectroscopy (SERRS), evanescent wave fluorescence (TIRF) and 
5 phosphorescence may be mentioned, as well as waveguide interferometers, 
waveguide lealdng mode spectroscopy, reflective interference spectroscopy 
(Rl^), transmission interferometry, holographic spectroscopy, and atohiic force 
microscopy (AFR). 

Commercially available biosensors include the BIACORE® system 
10 instruments, manufactured and marketed by Biacore AS, Uppsala, Sweden, 
which are based on surface plasmon resonance (SPR) and permit monitoring of 
surface binding interactions in real time between a bouiid ligand and an analyte 
of interest. 

While in the detailed description and Examples that follow, the present 
15 invention is illustrated In the context of SPR spectroscopy, and more partlculariy 
the BIACORE® system, it is to be understood that the present invention Is not 
limited to this detection method. Rather, any affinity-based detection method 
where an analyte binds to a ligand immobilised on a sensing surface may be 
employed, provided that a change at the sensing surface can be measured 
20 which is quantitatively indicative of binding of the analyte to the immobilised 
ligand thereon. 

• The phenomenon of SPR is well known, suffice it to say that SPR arises 
when light is reflected under certain conditions at the interface between two 
medfa of different refractive indices, and the interface is coated by a metal film, 

25 typically silver or gold. In the BIACORE® instalments, the media are the sample 
and the glass of a sensor chip which is contacted with the sample by a 
microfluidic flow- system. The metal film is a thin layer of gold on the chip 
■ surface. SPR causes a reduction in the intensity of the reflected light at a 
specific angle of reflection. This angle of niinimum reflected light intensity varies 

30 with the refractive index close to the surface on the side opposite from the 
reflected light, in the BIACORE® system the sample side. 
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A schematic illustration of the BIACORE® system is shown in Fig. 1 . 
Sensor chip 1 has a gold film 2 supporting capturing molecules 3, e.g. 
antibodies, exposed to a sample flow with'analytes 4 (e.g. an antigen) through a 
flow channd 5. Monochromatic p-polarised light 6 from a light source 7 (LED) is 
5 coupled by a prism 8 to the glass/metal interface 9 where the light Is totally 
reflected. The intensity of the reflected light beam 10 is detected by an optical 
detection unit (photodetector array) 11 . 

A detailed discussion of the technical aspects of the BIACORE® 
instrument and the phenomenon of SPR may be found in U.S. Patent No. 

10 5,31 3,264. iVIore detailed information on matrix coatings for biosensor sensing 
surfaces is given in, for example, U.S. Patent Nps. 5,242,828 and 5,436,161 . In . 
addition, a detailed discussion of the technical aspects of the biosensor chips 
used in connection with the BIACORE® instalment may be found in U.S. Patent 
No. 6,492,840. The full disclosures of the above-mentioned U.S. patents are 

1 s incorporated by reference herein. 

When molecules in the sample bind to the capturing molecules on the 
sensor chip surface, the concentration, and therefore the refractive ind®c at the 
• surface changes and an SPR response is detected. Plotting the response 
against time during the course of an interaction will provide a quantitative ' " • 

20 measure of the progress of the interaction. Such a plot is usually called a 
sensorgram. In the BIACORE® system, the SPR response values are 
expressed in resonance units (RU). One RU represents a change of 0.0001° In 
the angle of minimum reflected light intensity, which for most proteins and other 
biomolecules correspond to a change in concentration of about 1 pg/mm2 on 

25 the sensor surface. As sample containing an analyte contacts the sensor 
surface, the ligand bound to the sensor surface interacts with the analyte in a 
step refen-ed to as "association." This step is Indicated on the sensorgram by an 
increase in RU as the sample is initially brought into contact with the sensor 
surface. Conversely, "dissociation" nornially occurs when the sample flow Is 

30 replaced by, for example, a buffer flow. This step is indicated on the 
sensorgram by a drop in RU over time as analyte dissociates from the 
surface-bound ligand. 
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A representative sensorgram (binding curve) for a reversible interaction 
at the sensor chip surface is presented in Rg. 2, the sensing surface having an 
immobilised capturing mdecule, or llgand, for example an antibody, Interacting 
with a binding partner therefor, or analyte. in a sarnple. (The detection curves, 
or sensorgrams, produced by biosensor systems based on other detection 
principles mentioned above will have a similar appearance.) The y-axis 
indicates the response (here In resonance units, RU) and the x-axis Indicates 
the time (here In seconds). Initially, buffer Is passed over the sensing surface 
giving the baseline response A in the sensorgram. During sample injection, an 
increase in signal is observed due to binding of the analyte. This part B of the 
binding cun/e is usually referred to as the "association phase"; Eventually, a 
steady state condition is reached where the resonance signal plateaus at C (this 
state may, however, not always be achieved). At the end of sample Injection, ' 
the sample Is replaced with a continuous flow of buffer and a decrease in signal 
reflects the dissociation, or release, of analyte from the surface. This part D of 
the binding curve is usually referred to as the "dissociation phase". The analysis 
is ended by a regeneration step where a solution capable of removing bound 
analyte from the surface, while (Ideally) maintaining the activity of the llgand, is 
injected over the sensor surface. This is indicated in part E of the sensorgram. 
Injection of buffer restores the baseline A and the surface Is now ready for a ' 
new analysis. 

From the profiles of the association and dissociation phases B and D, 
repectlvely. Infomiation regarding the binding and dissociation kinetics is 
obtained, and the height of the resonance signal represents surface 
concentration (/.e., the response resulting from an interaction is related to the 
change In mass concentration on tiie surface). This will now be explained In 
more detail below. 

Assume a reversible reaction between an analyte A and a surface-bound 
(immobilised) capturing molecule, or llgand, B which is not diffusion or mass 
transfer limited and obeys pseudo first order kinetics: 
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A + B«*AB 

This interaction model (usually referred to as the Langmuir model), which 
assumes that the analyte (A) Is both monovalent and homogenous, that the 
llgand (B) is homogenous, and that all binding events are independent, is in fact 
s applicable in the vast majority of cases. 

The rate of change in surface concentration of analyte A (= rate of 
change in concentration of fornied complex AB) during analyte injection is the 
sum of the rates of the analyte A going on and off: 

^-Kum-kAAB] . (1) 

10 where [A] is the concentration of analyte A, [B] is the concentration of the llgand 
B, [AB] is the concentration of the reaction complex AB. ka is the association 
rate constant, and kd is the dissociation rate constant. 

After a time t, the concentration of unbound ligand B at the surface is 
[Bj] - [AB], where [Bj] is the total, or maximum, concentration of ligand B. 

15 Insertion into Equation (1) gives: 

^ = kM[Brhm-kAAB] (2) 

In tenns of detector response unite (AB is detected), this can be 
expr^sed as 

^''KCiR„^-Ii)-k,R (3) 

20 where R is the response at time t in resonance units (RU). C is the Initial, or 
bulk, concentration of free analyte (A) in solution, and R^ax "S the response (in 
RU) obtained if analyte (A) had bound to all ligand (B) on the surface. 
Rearrangement of Equation (3) gives: 

-^ = k„CR^-{k„C + k,)R (4) 

25 where R is the response in resonance unite (RU). In Integrated fonn, the 
equation is: 
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Now, according to equation (4), if dR/dt is plotted against the bound 
analyte concentration R, the slope is kaC +iC(j and the vertical intercept is 
' l^al^max^- If ^® ^^^^ concentration C Is known and Rmax t)een determined 
(e.g. by saturating the surface with a large excess of analyte), the association 
5 rate constant kg and the dissociation rate constant can be calculated. A 
more convenient metiiod is, however, fitting of the integrated function (5), or 
numerical calculation and fitting of the differential Equation (4), preferably by 
means of a computer program as will be described below.- 
The rate of dissociation can be expressed as: 

dR 

»0 -J = -M (6) 

and in integrated fonn: " 

R = H^e-'^ (7) 

wtiere Rq is tiie response at the beginning of tlie dissociation phase (when the 
huffei- wash of the surface starts). 
IS Equation (6) may be linearized: 

• ' 4f .. . . • 

and a plot of ln[R/Ro] versus t will produce a straight line with 0ie slope - -kd. 
More conveniently, however, the dissodatldn rate constant l<d is detennined by ' 
fitting the ^ponehtial rate equation (7). 
20 Affinity is expressed by the association constant Ka ~ ka/Rd, or the 

dissociation constant (also referred to as the equilibrium constant) Kq = kd/ka- 

The association constant Ka may altematively be determined from 
Equation (3), where dR/dt = 0 at equilibrium, giving: 

k.R^^WR^-R^) (9) 

25 Where Req is the detector response at equilibrium. Since ka/kd - Ka, insertion 
in Equation (9) and rearrangement gives: 

^ = -K,R^+K,R^ (10) 

If binding reactions are performed at multiple concentrations, the data 
may either be fitted or Req/C may be plotted versus Req which gives the slope 
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= -Ka- Such an equilibrium analysis may be perfomied when the rates of 
association and dissociation are too rapid to measure accurately. 

To obtain reliable kinetic constants, the above described analysis is 
usually repeated for a number of different analyte concentrations and, suitably, 
5 also at at least one other ligand density at the sensor surface. 

Software for the analysis of kinetic and other biosensor data is 
commercially available. Thus, for example, evaluation of kinetic data produced 
by the BIACORE® instruments is usually perfomied with the dedicated 
BIAevaluation software (supplied by Biacore AB, Uppsala, Sweden) using 

10 numerical integration to calculate the differential rate equations and non-linear 
regression to fit the kinetic parameters by finding values for the Variables that 
give the closest fit reducing the sum of squared residuals to a minimum. The 
residuals are the difference between the calculated and the experimental curve 
at each point, squared residuals being used to weight equally deviation's above 

15 and below the experimental curve. The sum of squared residual is expressed by 
Equation (11): 

^S-t(r,-rj' (11) 

where S is the sum of squared residuals, rf is -the fitted value at a given point, 
and Tx is the experimental value at the same point. 
20 - For example, for the molecular interaction described above, such 

software-assisted data analysis is performed by. after subtracting bacl<ground 
noises, making an attempt to fit the above-mentioned simple 1:1 Langmuir 
binding model as expressed by Equations (5) and (7) above to the 
measurement data. 

25 Usually the binding model is fitted simultaneously to multiple binding 

curves obtained with different analyte concentrations C (and/or with different 
levels of surface derivatization Rmax)- This is referred to as "global fitting", and 
based on the sensorgram data such global fitting establishes whether a single 
global kg or kj will provide a good fit to all the data. The results of the 

30 completed fit is presented to the operator graphically, displaying the fitted 
curves overlaid on the original sensorgram curves. The closeness of the fit is 
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also presented by the chi-squared (x^) value, a standard statistical measure. 
For a good fitting, the chi-squared value Is in the same magnitude as the noise 
in RU^. Optionally, "residual plots" are also provided which give a graphical 
Indication of how the experimental data deviate from the fitted curve showing 
5 the difference between the experimental and fitted data for each curve. The 
operator then decides if the fit is good enough. If not, the sensorgram or 
sensorgrams exhibiting the poorest fit are excluded and the fitting procedure Is 
run again with the reduced set of sensorgrams. This procedure is repeated until 
the fit is satisfactory. 

10 - Sometimes, the above-mentioned 1:1 binding reaction model will not be 
valid, which requires the data set to be reanalysed using one or more other 
reaction models. Such alternative models may include, for example, a one to 
one reaction influenced by mass transfer (mass transport), two parallel 
independent one to one reactions, two competing reactions, and a two state 

IS reaction. Parallel reactions can occur when the immobilised ligand is ' 

heterogeneous, whereas a heterogenous analyte may give rise to-competing 
reactions. A two state reaction indicates a confonnatlon change that gradually 
leads to a more stable complex between ligand and analyte. For differential rate 
• equations reflecting these alternative reaction models, It may be refened to, for 

20 example, Karlsson, R., and Fait, A. (1 997) J. Immunol. Methods 200, 121-1 33 
(the disclosure of which Is Incorporated by reference herein). For a more 
comprehensive description of curve fitting with regard to' the BIACORE® 
system, it may be refen-ed to the BIAevaluatlon Software Handbook (Biacore 
AS, Uppsala, Sweden) (the disclosure of which is Incorporated by reference 

25 herein). 

As regards mass transport, transport effects will influence the kinetics of 
binding when the reaction rate Is fast compared to the rate of transport. For a 
reaction where analyte (A) binds to immobilised ligand (B), binding with mass 
transport limitation may be represented by the reaction fomiula: 

2° Asolution <=> Asurface+ B <=> AB 

The differential equations describing the binding interaction will therefore 
include tenns for mass transfer of analyte to the surface. For a flow cell, a "two- 
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compartmenf model consisting of a set of coupled ordinary differential 
equations and described in inter alia l\/lyszka, D. G., et al. (1998) Bioptiys. J. 75, 
583-594, and Shanlc-Retzlaff, IVI. L. and Sligar, S. G., supra (the relevant 
disclosures of which are incorporated by reference herein) is considered to give 

5 a reasonable description of the binding l<inetics when the data are influenced by 
mass transport. In this model, the flow cell is assumed to be divided into two 
• compartments, one in which the concentration of analyte is constant; and a 
second near the sensor surface where the analyte concentration depends on 
the mass transport rate, the surface density of ligand. and the reaction rate 

10 constants. 

For the interaction of a monovalent analyte (A) reacting with an 
immobilised monovalent ligand (B), this model may be represented by the 
following two differential equations (12) and (13): 

dA 

Vi-^^S{-l^,A{Br-AB) + lc,AB^]c,{C-A)) . (12) 
at 

IS ^ = KA{Br-AB)-KAB (13) 

where Vj Is the volume of the inner compartment, S is the surface area of the 

flow chamber, is the transport coefficient describing diffusive movement of 

analyte between the compartments, Bj is the total ligand concentration, A is the 

bulk concentration of free analyte, C is the injection (I.e. initial) analyte ■ 
20 concentration, AB Is the concentration of complex AB (= surface density of 

bound analyte). and ka and k^ are the association and dissociation rate 

constants, respectively. 

It is understood that the use of this two-compartment model analyzing 

the kinetic globally will extend the range of reaction rates that can be 
2S determined with Instruments like the BIACORE® and has also been 

incorporated In the above-mentioned BIAevaluatlon analysis software (available 

from Biacore AB, Uppsala, Sweden). 

Theoretical interaction data for different Interaction models may be 

simulated using a simulator, such as e.g. the simulator which Is built-in in the 
30 above-mentioned BIAevaluatlon analysis software. In this way, the effects of 

varying experimental parameters may be explored, thereby helping to 
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determine which interaction conditions that may be crucial In distinguishing 
between alternative kinetic models. 

The invention 

s ' As mentioned above, determining kinetic parameters using, for example, 

a BIACORE® instrument, conventionally comprises monitoring the interaction of 
analyte with immobilised ligand for a plurality of different analyte concentrations 
with regeneration, or renewal, of the immobilised ligand before each 
measurement with a new analyte concentration; A kinetic model for the 

10 interaction is then fitted globally to the collected binding data, usually in the fomi 
of sensorgrams, to give the kinetic parameters. 

According to the present invention, such a flow cell-based analysis for 
- detemiining kinetic parameters may be substantially improved and speeded up 
by a titration type, or "sequential injection", procedure wherein the ligand- 

15 ' supporting surface is successively contacted, in one and the same analytical 
dycle, with the different analyte concentrations to produce a continuous 
sensorgram. A kinetic interaction model allowing for mass transport limititation, 
such as the two-compartment model mentioned above, is then fitted globally to 
the entire sensorgram to calculate the kinetic parameters. In addition to being 

20 considerably less time-consuming, this titration type method reduces the 
amount of experimental data to be evaluated and eliminates the risk of 
regeneration conditions destroying Immobilised ligand. ^ 

The order of introducing the different analyte concentrations is not critical 
to the success of the inventive procedure. Rather, the analyte concentrations 

25 may, in fact, be introduced in any order. For example, the analyte 

concentrations may be successively increased, or successively decreased, or 
may alternately be higher and lower, etc. Further, the same concentration of 
analyte may be introduced repeatedly, if desired. 

The injections of analyte, which preferably may be short, for a 

30 BIACORE® or similar instrument e.g. in the order of 30-60 seconds, may be 
interrupted when an increased analyte concentration only leads to a marginally 
increased response. This may advantageously be effected using so-called 
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adaptive software. The formed complex is then preferably allowed to dissociate 
for a longer period of time, e.g. from, about 5 to 60 minutes in the context of a 
BIACORE® or similar instalment, to add adequate dissodation data for the 
kinetic evaluation. Such a dissociation period may, optionally, instead (or 

5 additionally) be performed between any two analyte injections in the cycle. 

It is to be noted that the design of the above-mentioned BIACORE® 
instruments, for example, requires that each injection of analyte be followed by 
injection of running buffer, pemriitting partial dissociation of the bound analyte, 
before the next concentration of analyte is injected. It is understood that such 

10 alternating analyte-free liquid volumes, for example with about one minute 
intervals, will add dissociation data to the those obtained by the above- 
mentioned (preferably longer) dissociation period. Provided that sufficient 
dissociation data is obtained by these alternating injections of analyte-free fluid, 
the longer dissociation period may even be dispensed with. Intermediate buffer 

15 injections may, of course, also be used," even if a sensor instrument design as ' 
such permits sequential injections of analyte without buffer Injections between 
them. 

At least for some applications, it is prefen-ed that the flow is substantially 
constant during the entire experimental cycle. 

20 The sequential injection procedure may optionally be repeated at least 

once with the same or another sequence of analyte concentrations, and ail tiie 
collected binding data may then be fitted simultaneously to the kinetic 
Interaction model. 

If desired, binding data from one or more sequential experiments 

25 performed as as described above may be combined with binding data from one 
or more conventional type single concentration experiments, and the combined 
data are Uien fitted to the kinetic interaction model. Thus, generally, the total 
data set that is globally fitted to the kinetic interaction model may comprise any 
combination of binding curves including one or more "sequential" binding curves 

30 and one or more single concentration binding curves. 

To facilitate the understanding of th& invention, reference is now made to 
Fig. 3 which shows an exemplary sensorgram (as obtained with the above- 



18 



mentioned BIACORE® instrument) for a titration witli successively increased 
concentrations of an analyte (A), reacting with an Immobilised ligand (B) to form 
a complex (AB) which gives rise to the detected SPR response. This titration is 
described In more detail In Example 1 below. Referring to the sensorgram, 
5 analyte concentration 80 nM was injected between 404-464 seconds, 240 nM 
between 523-583 seconds, 730 nM between 642-702 seconds, and 2200 nM 
between 760-820 seconds. . 

Before analysing the binding curve, start and end times for each injection 
of analyte is identified, and the base line is adjusted to zero level. Generally, the 

10 binding curve may then be analysed using an interaction model based on time 
dependent analyte concentration and allowing for mass transport limjtation, 
such as the above-mentioned two-compartment model. If necessary, this model 
may be adapted to other Interaction mechanisms, such as e.g. conformation 
change (A + B AB o AB*). competing reactions (A1 +B A1B; A2 + B <-> 

15 A2B). parallel reactions (A + B1 ^ AB1 ; A +B2 <^AB2), or bivalent analyte (A + 
B AB; AB + B <-> AB2). Such models are described In more detail In, for 
example, the BIAevafuation Software Handbook mentioned above. 

Assuming, for example, 1:1 binding, the following differential equation 
may represent the time dependent concentration of analyte A capable of 

20 binding to immobilised ligand B. 
dA 

where A is the surface concentration of free analyte A, kt is a transport 
coefficient, C-| to Cp are the different bulk concentrations of analyte, B is the 
concentration of Immobilised ligand (ligand density), AB Is the concentration of 

25 formed complex (=bound analyte), and kg and kj are the association rate 
constant and dissociation rate constant, respectively. 

If at the beginning of the experiment the analyte concentration, A, at the 
surface is assumed to be zero, the concentration of ligand, B. is assumed to be 
Rmax» 9nd the concentration of formed complex, AB, is assumed to be zero, 

30 the following functions, here written in a form compatible with the above- 
mentioned BIAevaluation software, describe the interaction system. 
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AB+Rir$1+RI2*$2+RI3*$3+RI4*$4+RI5*$5: 
$1=(sign(t-tOn1)-(slgn{t-tOff1))/2: 
$2=(sign(t-tOn2)-(slgn(t-tOff2))/2; 
$3=(sign(t-tOn3)-(slgn(t-tOff3))/2; 
5 $4=(sign(t-tOn4)-(sign(t-tOff4))/2; 

$5=kt*($1*conc1+$2*conc2+$3*conc3+$4*conc4-A): 
$6=ka*A'B-kd*AB: 
A=$5-$6|0: 
B=-$6|Rmax; 
10 AB=$6|0 

In the above model, sign functions $1 to $4 (being 1 or 0) indicate when 
a certain analyte concentration Is Injected (tOn1 is the time when analyte 
concentration 1 is injected, tOff1 is the time when analyte concentration 1 1s 
stopped, etc). Functions $5 and $6 describe together the variation of the 
15 surface concentration of analyte over time, $6 describing how the ligand density 
and complex varies with time. The function on the first line con«sponds to the 
total change in signal from binding AB and bulk effects present during each 
Injection. By numerical Integration of these equations and global fitting to all 
data points of the whole sensorgram in Rg. 3, using e.g. BIAevaluation software 
20 (Blacore AB, Uppsala.'Sweden) adapted to the present invention, the kinetic 
rate constants ka and kd can be obtained. Mathematical modelis for the other 
interaction mechanisms may readily be devised by a person skilled in the art. 

The result of the fit is shown In f ig. 4, where the solid "non-rippling line" 
Indicates the adapted curve for best fit for the concentration range 80-2200 nM, 
25 kg and k^ being constrained to single values of 1.0*105 M-1s-1 and 6.6*10-4 s- 
1, respectively, as described in Example 1 below. 

With the above model, kinetic and affinity data can thus be obtained in 
spite of transport limitations. Since, as mentioned above, the method is • 
insensitive to the injection order of analyte, analysis of relatively low analyte 
30 concentrations may be performed when binding is close to saturation or above 
• the equilibrium response, thereby maximizing the kinetic infonmation in the 
sensorgram where transport effects are small. 
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Data from binding analyses at two or more different ligand densitites may 
advantageously be fitted simultaneously. 

Measurement of concentration simultaneously with kinetic rate constants 
may be perfomied, using the above titration procedure and data fitting, provided 

5 that the transport coefficient kt is known and kept constant. Preferably, binding 
data from at least two, and preferably more ligand densities are then used. 

A similar titration type approach as those above may also be used to 
determine (at least approximate) kinetic rate constants for a plurality of different 
analytes binding to the same ligand. In this case, different analytes are 

10 sequentially passed over the sensor surface rather than different analyte 

concentrations. The evaluation of the binding data is modified to pennit fitting of 
different rate constants to the resulting total binding curve. Such a procedure 
may, for example, be used to screen different binding partners towards the 
ligand to quickly find the best binder, e.g. drug candidates binding to a drug 

15 target, such as a receptor. Optionally, this approach may be combined with the 
use of different analyte concentrations. 

While it currently is preferred to use so-called label-free detection 
methods, Including those mentioned above, such as SPR. other detection 
techniques may. as already mentioned above, of course, also be contemplated, 

20 such as those based on the detection of a label, e.g. a radlolabel, fluorophore 
Onduding surface fluorescence), chromophore, chemllumlnescer, a marker for 
scattering light, an electrochemically active marker, a themnoactive mariner, etc. 

» 

The methods of the invention, including titration (I.e. sequential analyte 
injections) and evaluation with time dependent analyte concentration, are. as' 

25 already mentioned above, suitably reduced to practice In the form of a computer 
system running software which implements the different steps of Oie procedure. 
Preferably, the titration Is adaptive as mentioned above. The Invention also 
extends to computer programs, particularly computer programs on or In a 
canler, adapted for putting the titration procedure of the invention into practice. 

30 The cam'er may be any entity or device capable of carrying the program. For 
example, the canler may comprise a storage medium, such as a ROM, a CD 
ROM or a semiconductor ROM, or a magnetic recording medium, for example a 
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floppy disc or a hard disk. The carrier may also be a transmissible carrier, such 
as an electrical or optical signal which may be conveyed via electrical or optical 
cable or by radio or other means. Alternatively, the carrier may be an integrated 
circuit in which the program is embedded. 

5 Software like that outlined above may also be used to generate 

theoretical interaction curves for different interaction models similariy as 
described further above for'the conventional single injection analysis 
procedures. By generating data for a complex model, for instance a* 
conformation change model, varying parameters such as injection time and 

10 analyte concentration and then fitting this data to both a 1 :1 binding model and 
a confomriational change model the residuals can Identify conditions for optimal 
designs of kinetic experiments. 

In the following Examples, various aspects of the present invention are 
disclosed more specifically for purposes of illustration and not limitation. 

15 

EXAMPLES ' • 
Instrumentation 

A BIACORE® S51 instrument (Biacore AB. Uppsala. Sweden) was used. 
This Instrument has two Y-type flow cells which allow a dual flow of fluids over a . 

20 sensor chip surface, so-called hydrodynamic addressing, as described in, for 
example, EP-B1-1021703 (the disclosure of which is Incorporated by reference 
herein). The instrument uses three parallel detection spots on the sensor chip. 
As sensor chip was used Series S Sensor Chip CMS (Biacore AB, Uppsala, 
Sweden) which has a gold-coated surface with a covalently linked 

25 carboxymethyl-modified dextran polymer hydrogel. The output from the 

instrument is a "sensorgram" which is a plot of detector response (measured in 
"resonance units", RU) as a function of time. An increase of 1000 RU 
conresponds to an increase of mass on the sensor surface of approximately 1 
ng/mm2. 

30 

Example 1 

Binding of a kinase inhibitor to immobilised kinase 
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Anti-histidine antibody (Qiagen 34660; Qiagen, Venlo. Netherlands) was 
diluted to 20 Dg/mi in 10 mM acetate buffer pH 5.0. and 13000 RU of antibody 
was immobilised to a Series S Sensor Chip CMS using amine coupling 
according to the manufacturer's instructions. The flow was 10 ^l/min. IHistidine 
5 tagged kinase (proprietary) was then captured and cross-linked to the sensor 
chip surface by injecting 30 Dg/ml of the kinase in 10 mM PBS for three minutes 
followed by 2.5 min injections of EDO/NHS and ethanolamine to stabilise the 
kinase on the surface. With this procedure 4000 RU of kinase were 
immobilised. 

10 A kinase inhibitor, staurosporin, was injected at concentrations of 80, 

240, 730 and 2200 nIVl at a flow of 30 jil/min. The assay buffer was 50 mM Tris 
pH 7.5, 1 50 mM NaCI. 1.0 mM MgCl2 and 3% DMSO. Overlay plots of obsen^ed 
and fitted data for the injection sequence are demonstrated in Figure 3 
- (observed) and Figure 4 (observed plus fitted). The association rate constant 

15 was calculated to 1,0*10^ M"^ s"\ and the dissociation rale constant kd tq 
6.6*1 0"»s-^- 

Example 2 

Binding of a set of carbonic anhvdrase inhibitors to immobilised carbonic - 

20 anhvdrase 

Carbonic anhydrase (Sigma C3640), 20 bg/ml in 1 0 mM acetate buffer 
pH 5.0. was immobilised to a Series S Sensor Chip CM5 using amine coupling 
according to the manufacturer's instructions, (flow 10 nl/min) to obtain sensor 
surfaces with approximately 2000 and 700 RU, respectively, of immobilised 

25 carbonic anhydrase. Kinetic analysis was perfomied by injecting inhibitors in 10 
mM PBS buffer with 0.005% P20 and 5% DMSO at a flow of 30 ^il/min. The 
inhibitors were acetazoleamide (Sigma A601 1), azosulfamide (Sigma A2759), 
and benzenesulfonamide (Aldrich 10,814-6), Overlay plots demonstrating 
observed and fitted data for the two ligand densities are presented in Figure 5 

30 (acetazoleamide). Figure 6 (azosulfamide) and Figure 7 (benzenesulfonamide). 
In Figure 5, data represent two repeats on each ligand density, and ka was 
calculated to 1,8*10' M"^ s-\ and k^ to 4,0*1 0'^ s'\ In Figure 6. kg was 
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calculated to 2,9*10'* M"^ s'\ and to 8,1*10"^ s'\ and in Figure 7 kg was 
calculated to 1,6*10* M*^ s'\ and to 1,4*10"^ s'^- 

Example 3 

Binding of a high affinity Inhibitor to Immobilised thrombin 
Thrombin (Sigma T 1063). 20 Dg/ml In 10 mM acetate buffer pH 5.0, was 
Immobilised to a Series S Sensor Chip CMS using amine coupling according to 
the manufacturer's instructions (flow 10 jil/min), resulting In sensor surfaces 
with approximately 2300 and 600 RU, respectively, of immobilised thrombin. 
Kinetic analysis was performed by injecting a high affinity thrombin inhibitor, 
melagatran (a gift from Johanna Peinum, AstraZeneca, Mdlndal, Sweden), in 10 
mM PBS buffer with 0.005% P20, 5% DMSO and 3.4.mM EDTA at a flow of 30 
(il/mln. Overi^ay plots demonstrating observed and fitted data are presented in 
Figure 8. The calculated kg was 2.0*10' M'^ s'\ and kjj 1 .8*10"^ s"\ 

Examples 1 to 3 above demonstrate thpt a wide range of rate constants 
can be detemilned with the sequential Injection procedure. In these E)«imples. 
kg values range from 2,9*10* to 2.0*10'' M'^ s"^ and kd values range from 0.14 to 
6.6 10"* s'\ The Examples also illustrate that kinetic data can be obtaine<? using 
one or two llgand densities- and that the injection order can be mixed and not 
only in the order from low to high concentration. 

Example 4 

Combined sequential and single injections 
In Rgure 9, simulated data representing sequential and single Injections 
are combined in the evaluation where all binding curves are fitted 
simultaneously. The figure shows an overlay plot of four binding curves and 
fitted data. The three curves to the left represent sensorgrams with single 
Injections of analyte at 51200 nM but with varying injection times, 10 s, 30 s and 
500 s, whereas the curve to the right represents sequential injections of analyte 
from 800 to 12800 nM. This illustrates that sensorgrams with varying numbers 
of injections and injection times can be analysed. 
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it is to be understood that the invention is not limited to the particular 
emt)odiments of the invention described above, but the scope of the invention 
will be established by the appended dalms. 
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CLAIMS 



1 . A method of determining kinetic parameters for a reversibie 
molecular Interaction between a iigand immobilised to a soiid support surface 
and a binding partner to the iigand in solution, comprising the steps of: 

a) sequentially, v^rithout intemiediate regeneration or renewai of the 
immobiiised Iigand, flowing a plurality of fluid volumes containing different 
l^nown concentrations of the binding partner over the solid support surface to 
permit association of binding partner to the immobilised Iigand, 

b) flowing over the solid support surfoce a fluid volume free from binding 
partner to permit dissociation of binding partner from the Iigand, 

c) monitoring during steps a) and b) the momentary amount of binding 
partner bound to the soiid support surface related to time and solution 
concentration of binding partner and collecting the binding data, and 

d) detemilning the Icinetic parameters by globally fitting a predetermined 
kinetic model for the interaction between the binding partner and the 
immobiiised Iigand to the collected binding data, which model allows for mass 
transport limitation at the solid support surface. 

2. The method according to claim 1 , wherein steps a) to c) are 
repeated at least once with the same or another sequence of binding partner 
concentrations, and all the collected binding data are fitted simultaneously to 
the kinetic model in step d). 

3. The method according to claim 1 or 2, wherein sequential binding 
data obtained from steps a) to c) in claim 1 are combined with at least one set 
of binding data for a single concentration of binding partner, and the combined 
binding data are fitted In step d) to the kinetic model. 

4. The method according to claim 1 . 2 or 3. wherein step b) in claim 

1 is performed after flowing the last fluid volume containing binding partner over 
the suriiace. 
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5. The method according to any one of claims 1 to 4, wherein step b) 
in claim 1 is perfomned between two fluid volumes containing binding partner. 

6. The method according to claim 4 or 5, wherein each volume 
containing binding partner is followed by a fluid volume free from binding 
partner. 

7. The method according to claim 6, wherein one of the fluid volumes 
free from binding partner, preferably the last fluid volume, Is flowed over the 
solid support surface for a substantially longer time than the other fluid volumes 
free from binding partner. 

8. The method according to any one of claims 1 to 7. wherein the 
plurality of fluid volumes containing binding partner include fluid volumes with 

• the same concentration of binding partner. 

9. The method according to any one of claims 1 to 7, wherein the 
.fluid volumes containing binding partner contain successively Increasing 
concentrations of binding partner. 

1 0. The method according to any one of claims 1 to 9, wherein step a) 
in claim 1 is stopped when an increased concentration of binding partner results 
in a predetemiined marginal increase in the amount of bound binding partner at 
the sensor surface. 

1 1 . The method according to any one of claims 1 to 1 0, wherein the 
Idnetic model is a two-compartment model. 

1 2. The method according to any one of claims 1 to 1 1 , wherein the 
kinetic model Is adapted to an interaction mechanism selected from 1:1 binding, 
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parallel reactions, competing reactions, confonnation change and bivalent 
analyte. 

13. The nriethod according to any one of claims 1 to 12, wlierein step 
a) in claim 1 1s performed at at least one other density of llgand at the solid 
support surface. 

14. The method according to any one of claims 1 to 13, wherein tiie 
flow over tiie solid support surface is maintained essentially constant. 

15. A method of determining the concentrations of binding agent in a 
plurality of samples, wherein tiie binding agent Is capable of interacting with a 
llgand immobilised to a solid support surface, comprising the steps of: 

a) sequentially, wittiout intermediate regeneration or renewal of the 
immobilised ligand, flowing a plurality of different sample volumes containing ' 
binding agent over the solid support surface, ■ 

b) flowing over the solid suppoirt surface having binding agent bound 
thereto a fluid volume free from binding agent to pemiit dissociation of binding 
•agent from the ligand, 

c) monitoring during steps a) and b) tiie momentary amount of binding 
agent bound to the solid support surface related to time and solution 
concentration of binding agent and collecting the binding data, and 

d) detennlning tiie concenti-ations of binding agent In the different sample 
volumes by globally fitting a predetennined kinetic model for the interaction 
between tiie binding agent and tiie Immobilised llgand to tiie collected binding 
data, which model allows for mass taransport limitation at the solid support 
surface, and wherein the mass fransport coefficient is known and kept constant. 

16. The method according to claim 15, wherein also kinetic 
parameters are detennlned. 
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1 7. A method of determining kinetic parameters for reversible 
molecular Interactions between a ligand immobilised to a solid support surface 
and a plurality of different binding partners to the ligand in solution, comprising 
the steps of: 

a) sequentially, without intermediate regeneration or renewal of the 
immobilised ligand, flowing a plurality of fluid volumes containing different 
binding partners over the solid support surface to penmit association of binding 
partners to the immobilised ligand, 

b) after the addition of each fluid volume containing binding partner, 
flowing over the solid support surface a fluid volume free from binding partner to 
permit dissociation of binding partner from the ligand, 

c) monitoring duFing steps a) and b) the momentary amount of binding 
partner bound to the solid support surface related to time and solution 
concentration of binding partner and collecting the binding data, and 

d) determining the kinetic parameters by globally fitting a predetermined 
kinetic model for the interaction between each binding partner and the 
immobilised ligand to the collected binding data, which model allows for mass 
transport limitation at the solid support surface. 

18. The method according to any one of claims 1 to 17, wherein the 
solid support surface is a sensor surface. 

19. The method according to any one of claims 1 to 1 8. wherein the 
detection at the soiid support surface is based on evanescent wave sensing, 
preferably surface piasmon resonance (SPR). 

20. The method according to any one of claims 1 to 19, wherein the 
detection at the solid support surface is based on detection of a label. 

21 . The method according to any one of claims 1 to 20, wherein step 
d) in claims 1,15 and 17 is computer-implemented. 
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22. The method according to any one of claims 1 to 16, wherein steps 
a) to d) in ciaims 1,15 and 17 are computer-implemented. 

23. An analytical system for detecting molecular binding interactions, 
comprising: 

(i) a sensor device comprising at least one sensor surface, detection 
means for detecting molecular interactions at the at least one sensor surface, 
and means for producing detection data representing the progress of the 
interactions with time, and 

(ii) data processing means for performing the steps' as defined in any one 
of claims 1 to 20. 

24. A computer program comprising program code means for 
performing the method of any one of claims 1 to 20. when the program is run on 
a computer. 

25. A computer program product comprising program code means 
stored on a computer readable medium or carried on an electrical or optical 
signal for performing the method- of any one of claims 1 to 20 when the program 
is run on a computer. 

26. A computer system containing a program for perfomiing the 
method of any one of claims 1 to 20. 

27 Use of the computer program according to claim 24 or the 
computer program product according to daim 25 for simulating theoretical 
interaction data. 
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ABSTRACT 



A method of determining kinetic parameters for a reversible molecuiar 
interaction between a ligand immobiiised to a solid support surface and a 
binding partner to tiie ligand in solution, comprises sequentially, without 
intermediate regeneration or renewal of the immobilised ligand, flowing a 
plurality of fluid volumes containing different Icnown concentrations of the 
binding partner over the solid support surfece, monitoring the momentary 
amount of binding partner bound to the solid support surface related to time and 
solution concentration of binding partner and collecting the binding data, and 
detemnining the kinetic parameters by glot»IIy fitting a predetermined kinetic 
model for the interaction between tiie binding partner and the immobilised 
ligand to tiie collected binding data, which model allows for mass transport 
limitation at the solid support surface. An analytical system for carrying out the 
metiiod, a computer program, a computer program product and a computer 
system for performing tiie metiiod are also disclosed. 
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